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1. Introduction 

In order to obtain a deep understanding of the Quark-Gluon-Plasma medium and of its prop¬ 
erties, it is of central importance to study in which way the individual excitations of the system 
are modified with increasing temperature. Viewed globally, the hadron spectrum does not appear 
to change very much compared to the zero temperature situation until T ^ Tc where the transi¬ 
tion to the deconfined phase takes place. This observation is based on the success of the Hadron 
Resonance Gas Model (HRG) in describing equilibrium properties of the medium, such as quark 
number susceptibilities for T < Tc- These predictions are in very good agreement with Lattice 
calculations (see e.g. [1]). 

Nevertheless, detailed information about the individual degrees of freedom and of their dy¬ 
namics at finite temperature is lacking. This information can be extracted from thermal Euclidean 
correlators which can be calculated at the percent level in Lattice QCD. Here, we extend our lattice 
study of [2], where we derived and tested a modified pion dispersion relafion, Eq. (1.1), obfained 
from an effeclive fheory, which is valid for T <Tc 

tUp = u{T)^Jml+tp^. (1.1) 

The screening mass lUji is defined by fhe inverse correlafion lengfh and ft)p is fhe posifion of fhe 
pole on fhe ft)-axis in fhe refarded correlator Gr{(0,ip) wifh pion quanfum numbers. Under fhe 
assumpfion fhaf one is sufficienlly close fo fhe chiral limif and for small values of fhe momenfum, 
a single paramefer u{T) defermines fhe full dispersion relafion of soff pion quasiparficles af finife 
femperafure (for defails on fhe derivafion of Eq. (1.1) see [2]). In fhe chiral limif u{T) corresponds 
fo fhe group velocify of a massless excifafion Vg = dco/dp = n < 1 (in nafural unifs). Af T = 0, 
SO(4) Euclidean symmefry implies nij^ = (Oq and u = 1. 

In fhis proceedings article, we fesf fhe dispersion relafion Eq. (1.1) on a single fhermal fine 
laffice ensemble below Tc wifh high sfafisfics and fwo dynamical degenerafe lighf quarks [3]. In 
addifion, we revisif an old mefhod for inverting infegral equations, fhe Backus-Gilbert mefhod 
[4, 5], which, fo our knowledge, is used here for fhe firsf fime wifhin QCD in order fo obfain 
model independenf esfimators for specfral funcfions wifh inpuf from Euclidean correlators. This is 
particularly imporfanf since real-fime excifafions of fhe sysfem are encoded in specfral funcfions 
(e.g. via fhe hydrodynamic descripfion and Kubo formulas [6]). 


2. Lattice estimators of w at p = 0 


By exploiting fhe paramefric dominance of fhe pion in fime-dependenf axial-charge densify 
correlafors fogefher wifh Ward Idenfifies arising from fhe PCAC relafion we derived fwo indepen¬ 
denf laffice esfimafors for fhe paramefer u{T) [2] fhaf are valid af any femperafure below Tc (see 
also [7,8]): 


4m2 Gp(xo,r,0) 


ml Ga{xo,T,0) 

xo=Pl2_ 


^ _ flmn _ 

2G,4(j8/2,r,0)sinh(M/m;rj3/2) ■ 


( 2 . 1 ) 

( 2 . 2 ) 
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The relevant correlation functions that enter the calculation are 


5‘^‘’Ga{xo,T,p) = 

[dhe‘P-^{A/,{x)A^,iO)), 

(2.3) 

5^‘’Gp{xo,T,p) = 

j d\ e'P * 

(2.4) 


where a,b are flavor indices of the adjoint representation and jS = l/T. The quark mass is 
defined via the PC AC relation [9]. The screening pion decay constant as well as the screening 
mass are defined through the long distance behavior of Euclidean correlators along spatial 
directions 


= jdxodW{A'i{x)Al{<d)) (2.5) 

5‘'’’G^p{x3,T) = j dxod^X^(F\x)P\<d)^ ^ 2 . 6 ) 

Notice that the parameter u{T) is a RG invariant quantity (so are and Uf) and therefore renor¬ 
malization constants cancel out. 


3. Lattice setup 

Measurements where performed on a 24 x 64^ lattice with two dynamical light quarks with a 
mass of niq = = 2GeV) = 12.8(l)MeV. The temperature is T = 1 jlAa = 169MeV and the 

spatial extent amounts to L = 64a = 3.1fm. In the chiral limit the critical temperature for Nf = 2 
is rc(0) Ri IVOMeV [10] and since Tc{niq) grows with niq our thermal ensemble lies below the 
phase transition. In addition, an effective zero-temperature 128 x 64^ ensemble with the same bare 
parameters is available to us through the CLS effort (labelled as 07 in [11]). Therefore we are able 
to compare in a straightforward manner our finite temperature results with the vacuum situation. 
Results are shown in Table 1 and a comparison with previous work is shown in Fig. 1. 


mn/T 

1.79(2) 



fn/T 

0.46(1) 

(Oo/T 

1.579(12) 

Uf 

0.76(1) 

(Oi/T 

2.88(3) 

Ujn 

0.74(1) 

fjT 

0.599(8) 

(Oo/T 

1.32(2) 

M2(r~o) 

1.01(6) 

fjT 

0.62(1) 




Table 1: Left: Summary of the results for the = 24 thermal ensemble. The value of cuo is calculated 
using Oo = Uminjc- In the same way = fnlum- Right: Summary for the vacuum ensemble 07. The values 
of ©0 and cui are extracted from the long range behavior of Eq. (2.1) and Eq. (2.2) with p = {Q,Q,2nn/L) 
and n = 0,1 respectively. The value of 1^(7 ~ 0) was extracted via a linear fit to the square of Eq. (1.1) 
with mji replaced by (Uo- All renormalization factors are included and the errors are purely statistical. All 
dimensionfull results are normalized with T = I /24a. 
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Figure 1: Comparison of the results obtained in [2] with this work [3]. Notice that in [2] the temperature 
scan was performed at a constant renormalized quark mass of — 2GeV) ^ 15MeV. 

4. Axial-charge density correlator at p 7 ^ 0 

The relation between the Euclidean correlator Ga{xq, r,p) and its spectral function is 

cosh(a)(j8/2-;co)) 

UUJ fj Vtu,p 
/O 


Ga{xo,T,p)= d(op^{co,p) 


(4.1) 


sinh(ft)j3/2) 

At finite temperature, the analysis of the correlator Ga{xo, r,p) is more involved than at zero tem¬ 
perature: only at sufficiently small quark masses and momenta, and not too small xq, is the cor¬ 
relator parametrically dominated by the pion pole. Therefore, we proceed by formulating a fit 
ansatz for the spectral function to take into account the non-pion contributions with the property 
that p'^ (ft) —^ ooj r, p) = 6 (co^ — 4m^ — p^) 2 ^ (p^ + 6 m^) (see e.g. [12]). By integration with the 
kernel spectral function and the correlator fit ansatze read 

\-e-^fe{(0-c), (4.2) 


Nr 


P («,P) = 4i(p)sinh(ftjj8/2)5(ftJ-ftJp)-hA2(p)^^ 


Nr 


-CXq 


Ga{xo,T,p) = Ai(p)cosh(ftjp(j8/2-xo))-hA2(p)^^ 


xo 


j8 -Xo 


(4.3) 


The fit involves only 4 parameters but leaving cUp free led to poorly constrained fits. Therefore, we 
fix fhe value of ft)p fo fhe predicfion of Eq. (1.1) wifh u = Um and p = {0,0,2jin/L) fo see whefher 
fhe dafa can be described in fhis way. Resulfs are shown in Table 2. In view of fhe ;^^-values, fhe 
dafa is consisfenf wifh fhis scenario. Moreover, fhe chiral effeclive fheory makes a predicfion for 
fhe residue of fhe refarded correlafor at co = (Op reading (see App. B of [3]) 


P^{(0,p)=f^{ml+p^)5{(0^-(ol) + 


Res(a)p) 


(4.4) 


which can be wriffen in terms of A\ (p) and a new parameter Z7(p) which parametrizes the deviation 
with respect to the chiral prediction 


Res(ft)p) = 2Ai(p)ci)pSinh(ftjpjS/2) 

= fn{ml + p^){l+b{p)). 


(4.5) 

(4.6) 
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In addition, we rescale the parameter A 2 (p) —)• A 2 = A 2 (p)/p^ (neglecting the quark mass which is 
parametrically small) whose natural value is 1. It turns out that indeed the parameter Z>(p) is small 
for n = 1 and the values of A 2 are of order one adding confidence to our description. 


n 

Ai/r3 


A 2 

c/T 

Res(mpJ/r4 

b 

Z^/d.o.f 

1 

2.95(4) X 10-' 

2.19(3) 

1.78(8) 

6.7(3) 

1.72(6) 

-0.08(3) 

0.06 

2 

1.40(5) X 10-' 

3.73(6) 

1.26(2) 

6.1(1) 

3.3(2) 

-0.39(4) 

0.15 

3 

4.9(3) X 10-2 

5.40(9) 

1.19(1) 

7.7(1) 

3.9(5) 

-0.65(4) 

0.35 

4 

1.7(2) X 10-2 

7.1(1) 

1.15(1) 

9.67(9) 

4.21(7) 

-0.78(3) 

0.49 

5 

4(1) X 10-2 

8.8(1) 

1.12(1) 

11.7(1) 

3(1) 

-0.89(3) 

1.04 


Table 2: Results of fits to the axial-charge density correlator at non-vanishing momentum p„. All errors 
quoted are statistical, and all renormalization factors are included. The quantity CUp/r is not a fit parameter; 
rather it is set to the value predicted by Eq. (1.1) with u{T) = u„, = 0.74(1). 


5. Backus-Gilbert method on Ga {xq, T, p) 


The Backus-Gilberl melhod (see [13,14] for more recenl sludies) is a lechnique for inverting 
integral equations like Eq. (4.1) in order lo direclly oblain information on p'^(ft),p). The idea is lo 
define an eslimalor 

p(S.p)=/Gco5(S.»)(;i|g^) (5.1) 

where 5{(d, (o) is called fhe resolufion funclion which can be wriflen in ferms of a sel of priori un¬ 
known coefficienls qi{6)) and fhe kernel functions Ki{(o) = ^ smoolh function 

peaked around (b which "filters" fhe Irue speclral funclion. The coefficienls qi{(b) are delermined 
by minimizing ifs widlh subjecf lo fhe condition lhal fhe area is normalized lo 1. During fhe process 
fhe malrix W need lo be inverted. 


Wij{(b) = X 


d(oKi{(o){(0 - (o)^Kj{(o) 


+ {\-X)Cov{Gi,Gj), 0<A<1 


(5.2) 


where Cov{Gi,G j) denotes Ihe /y-elemenl of fhe covariance malrix of Ga{xq.,T ,tp). The parameter 
X conlrols fhe Irade off belween slabilily and resolution power. For values of X close lo 1, fhe 
malrix is almost singular and the error on p{(b) rapidly increases. Reducing the value of X improves 
stability at the cost of deteriorating the frequency resolution. Finally the solution reads 


p(m,p)=£^,-(«)GA(4,r,p). (5.3) 

i 

Notice that p(d),p) is a completely model independent estimator in the sense that no ansatz was 
needed for its calculation. Fig. 2 shows the estimators for different values of p and one already sees 
the good agreement in the high-frequency region with tree-level predictions. We can further check 
the agreement of Res(ft)p) with the chiral prediction by making use of p((i),p). Assuming a 5-type 
excitation in we can write 

Res(«^,co)Bo = (54) 

5(0),cUp) 
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Figure 2: Estimators p{(0, p„)/7’^ for n — 0,1, 2,3 together with its error shown as a band (A = 0.002). The 
vertical colored dashed lines correspond to the locations of the expected positions of the poles Wp according 
to Eq. (1.1) with u{T) = The black horizontal lines correspond to the treelevel asymptotic values of 
p^{(0,p). All renormalization constants have been taken into account. Dimensionful quantities have been 
made dimensionless by the appropriate power of T = 1 /24a. 

where we use ft)p as input according to Eq. (1.1) with u{T) = u,„. This defines a function of co. The 
natural choice where Res(tUp, tu)BG is expected to be the best estimator for the true residue is at 
ft) Ri ft)p. By looking at Fig. 3 we see that approximately around this value the curve intercepts the 
grey band which represents the chiral prediction and therefore confirms fhe validify of Eq. (1.1) af 
leasf up fo p = 400MeV {n = 1). 



ui/T u/T 


Figure 3: The effective residue Res(ft)p, ft))BG dehned in Eq. (5.4). Left: No momentum induced, p = 0. 
Right: One unit of momentum induced, pi = {0,0,2k/L). The grey band is the expectation in terms of 
screening quantities. All renormalization factors are included. The errors arise from the statistical uncer¬ 
tainty. The values of ft)p are indicated by dashed vertical lines. 

6. Conclusions 

We have successfully tested the dispersion relation of Eq. (1.1) for the pion quasiparticle by 
performing direct fits to Ga/xq, T,p) (alternatively applying the Backus-Gilbert method) using the 
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value of the parameter u determined at zero momentum. Table 1 indicates that the pion mass 
"splits" at finite temperature into a lower pole mass and a higher screening mass. Our findings 
are in qualifative agreement with two-loop calculation in ChPT (see eg. [15,16]). Both suggest a 
violation of boost invariance due to the presence of the medium. Although = 24 is at the current 
state of the art concerning thermal ensembles, finite volume effects as well as lattice artifacts should 
be investigated. Our plans for the future are to simulate at the physical value of the light quark mass 
and/or to include the strange quark. 
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